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We investigate the phenomenon of quantum radiation - 
i.e. the conversion of (virtual) quantum fluctuations into (real) 
particles induced by dynamical external conditions - for an 
initial thermal equilibrium state. For a resonantly vibrating 
cavity a rather strong enhancement of the number of gener- 
ated particles (the dynamical Casimir effect) at finite temper- 
atures is observed. Furthermore we derive the temperature 
corrections to the energy radiated by a single moving mir- 
ror and an oscillating bubble within a dielectric medium as 
well as the number of created particles within the Friedmann- 
Robertson- Walker universe. Possible implications and the rel- 
evance for experimental tests are addressed. 

PACS: 42.50.Lc, 03.70.+k, ll.10.Ef, ll.10.Wx. 



I. INTRODUCTION 

Motivated by previously obtained results Q we gener- 
alize the canonical formalism adopted there towards the 
application to further scenarios and establish the deriva- 
tions in more detail: 

One of the main consequences of quantum theory is 
the existence of a non-trivial vacuum state. In contrast 
to the classical theory the quantum fields undergo fluc- 
tuations even in their state of lowest energy (the ground 
state) - the so-called vacuum fluctuations. These fluctu- 
ations have measurable consequences: E.g., if the fields 
are constraint by the presence of external conditions, the 
energy associated to these fluctuations (the zero-point 
energy) may change owing to the imposed external con- 
ditions. As a result the quantum field may exert a force 
onto the external conditions in order to minimize its en- 
ergy. The most prominent example for such a force is 
the Casimir [Q effect which predicts the attraction of two 
parallel perfectly conducting and neutral plates (i.e. mir- 
rors) placed in the vacuum of the electromagnetic field. 
The prediction of this striking effect has been verified 
experimentally with relatively high accuracy |3|,|| . 

A different - not less interesting - effect has not yet 
been rigorously verified in an experiment: The impact 
of the external conditions may also induce a conversion 
of the virtual quantum fluctuations of the field into real 
particles - the phenomenon of quantum radiation. As 
examples for such external conditions giving rise to the 
creation of particles we may consider moving mirrors, 
time-dependent dielectrics or gravitational fields. 

Various investigations have been devoted to this topic 
during the last decades, in order to mention only some of 
the most important initial papers in chronological order: 



In 1970 Moore || presented the first explicit calculation 
of the quantum radiation on the basis of two 1+1 di- 
mensional moving mirrors. In this pioneering work he 
exploited the conformal invariance of the scalar field in 
1+1 dimensions. Based on this result Fulling and Davies 
H presented a calculation of the radiation of a single 
moving mirror (again in 1+1 dimensions) and pointed 
out the close analogy to the Hawking radiation. In 1982 
Ford and Vilenkin [Q succeeded to develop a method for 
the calculation of the radiation generated by a moving 
mirror in higher dimensions, i.e., without exploiting the 
conformal invariance. 

For the experimental verification of the phenomenon of 
quantum radiation the scenario of a closed cavity might 
be most promising since one may exploit a resonance en- 
hancement in this situation. The particle production 
inside a resonantly vibrating cavity has already been 
considered by several authors, see e.g. [p|-p"7| as well as 
|ij|-f20| for reviews. 

However, most of the investigations of quantum ra- 
diation are restricted to the vacuum state, i.e. to zero 
temperature. But in view of an experimental verifica- 
tion it is essential to study the finite temperature effects. 
Realistic calculations of thermal effects on quantum ra- 
diation within the framework of quantum field theory of 
time-dependent systems at finite temperature are not yet 
available. 

The remedy of this deficiency is the main intention of 
the present article: In Sec. || we set up the basic formal- 
ism for the quantum treatment of external conditions at 
finite temperatures. The developed methods are applied 
in Section III to the scenario of a trembling cavity. In 
Sec. |y] we focus on the resonance case and derive the 
number of created particles. Another scenario giving rise 
to quantum radiation - a dynamical dielectric medium 
- is considered in Section In Sec. VI we demonstrate 
the flexibility of the canonical approach presented in Sec- 
tion H by calculating the finite temperature corrections 
to the particle production in yet another example sce- 
nario - the Friedmann-Robertson- Walker universe. We 
shall close with a summary, some conclusions, a discus- 
sion, and an outlook. 

Throughout this article natural units with 



% = c = G N — k B = £o — Mo — 1 



(1) 



will be used. The signature of the Minkowski metric is 
chosen according to g^ v = diag(+l, —1, —1, —1). 



1 



II. GENERAL FORMALISM 

The objective is to investigate quantized bosonic fields 
obeying linear equations of motion under the influence 
of external conditions. At asymptotic times \t\ f oo 
the external conditions, which are treated classically 
(not quantized), are assumed to approach a static Hq- 
configuration, where the asymptotic Hamiltonian Hq can 
be diagonalized via a suitable particle definition. Ini- 
tially the state of the quantum system is supposed to 
be described by thermal equilibrium at a given temper- 
ature T, which might be realized through the coupling 
to a corresponding heat bath. However, the coupling to 
the heat bath has to be switched off before the external 
conditions undergo dynamical changes in order to avoid 
relaxation processes (closed system) . In general the time- 
dependence of the external conditions cause the state of 
the quantum system to leave the initial thermal equi- 
librium. Accordingly, the calculation of the expectation 
value of relevant observables, e.g. the number of particles, 
before and after the dynamics may deviate. These differ- 
ences can be interpreted as particles that are created or 
even annihilated by the dynamical external conditions. 



A. Interaction picture 

The following calculations are most suitably performed 
in the interaction representation. Accordingly, the dy- 
namics of all operators X corresponding to observables 
are governed by the undisturbed Hamiltonian Hq 



dX_ 
~~dt 



H ,X 




(2) 



This Hamiltonian Ho describes the complete dynamics of 
the system at asymptotic times and can be diagonalized 
via a suitable particle definition 



H(\t\ t oo) = Ho 



E Q = ui Ni + Eo 



(3) 



where Eo denotes the (divergent) zero-point energy. In 
most of the following formulas we make use of a gener- 
alized sum convention and drop the summation signs by 
declaring that one has to sum over all indices that do 
not occur at both sides of the equation. Equations with 
the same index appearing at both sides are valid for all 
possible values of this index. 

The index / contains a complete set of quantum num- 
bers labeling the different particle modes, e.g. / = {k} 
or / = {cu,£,m} etc. The particle energies are given by 
Lj[. For a thermodynamical consideration we have to de- 
scribe the state of the field by the statistical operator 
g. In the interaction picture the time evolution of this 
density matrix is given by the von Neumann equation 



dp 



Hi,g 



(4) 



The perturbation Hamiltonian Hi governs the influence 
of the variation of the external conditions upon the quan- 
tized field. Note, that this equation describes the time 
evolution of a closed quantum system, i.e. no measure- 
ments (etc.) take place during the dynamics. It leads 
to an unitary time evolution operator and therefore does 
also not contain relaxation processes, etc. Measurements 
and relaxations would change the probabilities wa of the 
statistical operator 



g(t) = w A \9 A (t)) 



(5) 



and can be incorporated into Eq. (Q) via an explicit time 
derivative (dg/dt) cxp . 



B. Entropy 

Without an explicit time-dependence (dg/dt) cxp 
Eq. (Q) generates a unitary time evolution and hence, 
the microscopic entropy remains constant in time 



S = — Tr{ghig} = const. 



(6) 



Note, that a constant microscopic entropy arises also in 
classical mechanics where the time evolution is governed 
by the Liouville equation. By virtue of Liouville's the- 
orem the total time derivative of the phase space den- 
sity g vanishes and thus the classical microscopic en- 
tropy J dTghig remains constant as well. But introduc- 
ing the Boltzmann equation via averaging over multi- 
particle correlations it is possible to define an effective 
entropy which increases in general (_ff-theorem). 

An analogous procedure can be performed in quan- 
tum theory: In practice, a complete knowledge about 
a given quantum system can never be achieved. For- 
mally, this restriction defines a so-called observation level 
Q = {X} as a set of possibly relevant observables X (see 
H), where an averaging over all unknown and possi- 
bly irrelevant observables is understood. With respect 
to a given observation level Q one can introduce an ef- 
fective statistical operator §{gy such that it yields the 
correct expectation values (X) = Ti{grg\X} = Tr{gX} 

for all operators X G Q. The effective statistical opera- 
tor g^gy averages over all irrelevant observables X £ Q in 
order to maximize the effective entropy which is defined 
as S{gy = — Trjgjg;} In g\gy\- This effective entropy S^gy 
referring to a given observation level in general increases 
with time. 

Introducing the internal energy E = (: Ho :) as ob- 
servation level the corresponding effective entropy Se 
may increase under the influence of the dynamical ex- 
ternal conditions reflecting the fact that particles have 
been created. The physical meaning of Se respectively 
its change ASe may become evident, if one assumes some 
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energy-conserving relaxation process, e.g. mediated via a 
measurement after the dynamics has taken place, which 
thermalizes the system again at some higher, in principle 
measurable temperature Te = T + AT. For a photon 
gas we find in the limit of high temperatures respectively 
of large volumes 15 the following relations between the 
energy E, the effective entropy Se and the effective tem- 
perature T E : E = l3T E n 2 /15 and S E = 15 T E An 4 / 45 . 
For small disturbances the relative increase of the effec- 
tive temperature, the energy and the effective entropy 
behaves as AE/(4E) « AT E /T E » AS E /{ZS E ). 



C. Time evolution 

Within our approach the influence of the dynamics of 
the external conditions is represented by the perturba- 
tion Hamiltonian H\{t) governing the time-evolution of 
the statistical operator in Eq. (Q) . By means of the time- 
ordering operator T this equation can be integrated for- 
mally 



g(t T oo) = U g(ti -oo) £7* 



= r 



exp ( -i / dt' H x (t') 



g(t i -oo) 



exp [i / dt" Hx{t" 



The chronological operator T acting on two bosonic and 
self-adjoint operators X(i) and Y(t') is defined by 



T 



X(t)Y(t') =X(t)Y(t')Q(t-t') 

+Y(t')X(t)Q(t' - 1) , 



(8) 



and so on for more operators. Due to the Hermitian 
conjugation of the unitary time evolution operator U in 
Eq. (0), for which the position of all operators changes, it 
is convenient to introduce the anti-chronological operator 
Tt (cf. H) as well 



rt 



X(t)Y(t') 



T 



X(t)Y(t') 
Y(t')X(t)Q(t-t') 

+x(t)Y{t')e(t' - 1) 



(9) 



Combining both equations one obtains for two bosonic 
and self-adjoint operators 



|x(t),Y(i')} = T [x(t)Y(t')\ +T+ [x(t)Y(t') 



Y(t')X(t)+X(t)Y{t') 



(10) 



As we shall see below this property simplifies the calcu- 
lation of the quadratic response of the number operator. 



D. Canonical ensemble 

As stated in Section || we assume the quantum field to 
be initially at thermal equilibrium corresponding to some 
non- vanishing temperature T = 1/(3 > 0. For reasons of 
simplicity we restrict our further consideration to parti- 
cles, that do not exhibit another conserved quantity than 
their energy. 1 This assumption is correct for photons, 
but not for charged pions, for instance. We consider only 
bosons. For that reason the chemical potential vanishes 
and the energy E is the only one observable that has 
a fixed expectation value. Minimizing the microscopic 
entropy with this constraint generates the canonical en- 
semble 



£?o = g(t i -oo) = 



exp ~(3H 



Tr {exp (-Affo) } 



(11) 



There are several examples in quantum field theory where 
the canonical ensemble is not capable for describing the 
thermal equilibrium correctly, usually in combination 
with infinite volumes. In order to treat also such cases, 
we assume the system to be confined into a finite volume 
where the canonical ensemble applies, calculate the ex- 
pectation values (i.e. the trace), and consider the infinite 
volume limit afterwards. 



E. Response theory 



In general, the final expectation value of an observable 



X 



Tr{l£(iToo)} 
Tr{l U g(t i -oo) #t\ 



(12) 



cannot be calculated explicitly for non-trivial interaction 
terms Hi owing to the complicated structure of the cor- 
responding time-evolution operator U. For that purpose 
one has usually to introduce some approximations. One 
possibility is given by the perturbation expansion with 
respect to powers of the disturbance Hi . Assuming the 
perturbation Hamiltonian Hi to be small it is possible to 
expand the above expression in powers of Hi . Neglecting 
all terms of third and higher order in Hi one obtains the 
quadratic response 



X ) = Tr 



{Xgo} 



TriX 



g ,i / dtHi(t) 



+Tr\x / dtHi(t)g / dt' Hi{t') 



1 Otherwise, one may start with a grand canonical ensemble. 
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1. 



Tri XT 



dtH^t) / dt' Hi(t') 



L>() 



--Tr {Xg T^ 



dtH^t) / dt'Hi(t') 



Hi 



(13) 



Focusing on the investigation of the particle production 
the relevant observable is the number operator X = Nj. 
Due to [iVj, go] = the linear response vanishes and with 
the aid of Eq. ( jToj ) the quadratic response simplifies to 



N T ) = Tr 



{NiQo} 



-Tr ^ iVj 



dtHi(t), g 



dtH x {t) 



-O 



N, 



ANj + O 



Hi 



(14) 



The first term (Ni)o of the above expression denotes the 
initial particle content in the canonical ensemble which 
is given by the Bose-Einstein distribution. The second 
term ANi describes the particle creation or annihilation 
due to the presence of dynamical external conditions and 
will be calculated in the following. 



F. Particle production 



Note, that there is still a summation over the index J. 
One observes that merely the 5-term governs the pro- 
duction of particles and contains the vacuum contribu- 
tion (first addend). The U-term does not increase the 
total number of particles since it has the same structure 
as a classical master-equation (e.g. used for the deriva- 
tion of the if -theorem) , but it transforms particles from 
one mode into another and thereby also increases the 
energy. Investigating the high-temperature expansion 
of the Bose-Einstein distribution entering the equation 
above 



N, 



1 



o exp(/3cjj) — 1 



1 



(17) 



one observes that the temperature-independent contri- 
butions (the term —1/2) cancel. The same occurs in 
the static Casimir effect |23|-^f| and may be interpreted 
as a consequence of the scale invariance in the classi- 
cal limit. In the high-temperature limit the expectation 
value (Ni)o is linear in T. But special care is required for 
evaluating the number of created particles in Eq. ([II]) due 
to the remaining mode summation. Since the expansion 
in Eq. ( fL7|) has a finite radius of convergence its insertion 
into Eq. ( |l6| ) may cause some problems in performing 
the mode summation. Therefore the resulting number of 
produced particles or the total radiated energy may pos- 
sibly display another behavior a s show n i n Eq . (|17|). As 



we shall see later in the Sections III C and V A , the radi 



ated energy may be proportional to even higher powers 
in T. 



For quasi-free quantum fields obeying linear equations 
of motion the perturbation Hamiltonian can be expressed 
- up to an irrelevant constant - as a bilinear form of the 
fields. Expanding the fields into the time-independent 
creation/annihilation operators representing particles in 
the asymptotic regions the disturbance H\ can be cast 
into the rather general form 



1 



dtHi{t) = - [SjKd\cJ K +S*j K aja,K 

+UjKO { jdK + C . 



(15) 



The introduced matrices have to fulfill the conditions 
Sjk = Skj and Ujk = W^j because H\ is self-adjoint. 
The 5-term could be interpreted as a generator for a 
multi-mode squeezing operator and the W-term as a hop- 
ping operator. Now the matrices S and IA contain all in- 
formation about the dynamical external conditions that 
are relevant for the quadratic response. Inserting the 
general form of the disturbance in Eq. (15) into Eq. (JL4) 



after evaluating the traces the quadratic response of the 
number operator takes the form 



AJV> = \Sij\(l + {Nj) +(JV> 



\Uu\ 2 {( 



N 



o \ / o 



(16) 



G. Correlations 

The expectation value (iVj) of the number operator 
after the dynamical period does in general not display 
the Bose-Einstein distribution because the field is no 
longer in the thermodynamic equilibrium. But for cer- 
tain dynamics this expectation value could still behave 
as a thermal distribution corresponding to some effective 
temperature T c g. If particles are created during the dy- 
namics this effective temperature will be larger than the 
initial temperature T. Such a phenomenon may occur 
even for a vanishing initial temperature T = 0. This ef- 
fect can be explained within the thermo-field |2(|^7j for- 
malism: Measuring only single-particle observables does 
not reveal the complete information about the quantum 
system under consideration. Formally, this restriction 
defines an observation level Q pl| including all these 
single-particle observables (see the remarks in the pre- 
vious Section). The effective density matrix Q{g} may 
indeed be equal to the statistical operator of a canon- 
ical ensemble corresponding to some effective tempera- 
ture T c ff. However, the real density matrix g cannot be 
equal to this effective statistical operator g^g} because 
the microscopic entropy S — — Trjplnp} is conserved 
during a unitary time-evolution while the effective en- 
tropy S{g} = — Tr{^{0} In Q{gy} has been increased for 
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T < T c ff. Within the investigation of only single-particle 
observables one can never distinguish between the two 
statistical operators g and g{gy and therefore one can 
never find out whether the measured temperature repre- 
sents a real one (T) or an effective one (Toff)- For this 
purpose it is necessary to consider many-particle observ- 
ables. One suitable candidate is given by the two-particle 
correlation defined as 



C. 



JK 



NjN K 



N.j) N 



for J^K. (18) 



This quantity is particularly appropriate since it vanishes 
in the thermodynamic equilibrium 



NjN K 



for J ^ K . 



(19) 



The quadratic response of the correlation function can 
be evaluated as follows 



Cjk — 



dtHx(t) 



NjN K , / dtHj.it) 



Nj) AN K 
o 



N K ) ANj + O 



Hi . (20) 



As a simple example we may consider the vacuum case 
with T = 



Cjk — 



dt Hi(t)NjN K / dtHx{t)\Q) 



Hi 



(21) 



where the correlation is positive (at least in lowest order 
in Hi). 

But at finite temperatures the correlation may also as- 
sume negative values: E.g., in the case of a completely 
diagonal perturbation Hamiltonian we arrive at 



C 



JK 



Nj) q ANk-[Nk) o ANj 



O 



Hi 



(22) 



H. Bogoliubov transformation 

It might be illuminating to consider the relation of the 
previous investigations to the formalism based on the Bo- 
goliubov transformation. The initial and final creation 
and annihilation operators are connected through the Bo- 
goliubov coefficients via 



U^ajU = ajK&K + PjkoX 



K ■ 



(23) 



Switching from the interaction picture to the Heisenberg 
representation the expectation value of the number op- 
erator 



Tr 



[gNj] =Tr [g U^NjU] 

= (frNjU^ (24) 



can be expressed in terms of the Bogoliubov coefficients 



•Vy ) = \(3jk\ 2 + (Nk) (k/if I 2 + l/W) ■ (25) 



If we assume a completely diagonal Hamiltonian H — 
H n + Hi the associated time-evolution operator U fac- 
torizes 



(26) 



Accordingly the Bogoliubov coefficients simplify to 
oljk = ajSjK and (3 JK = PjSjk, respectively. Utiliz- 
ing the unitary relation, which assumes for the diagonal 



coefficients the simple form |cc/| 2 
at 



1 



1/3/ 



we arrive 



Nr ) = (N, 



\Pi\ 2 



1 + 2(N I 



(27) 



Hence - for a diagonal Hamiltonian - the number of cre- 
ated particles at finite temperature is simply given by the 
corresponding vacuum expression times a thermal factor 



AN 1 



AN, 



T=0 



1 



exp(/3w) 



1 



(28) 



It should be mentioned here that this result is not re- 
stricted to a particular order perturbation theory - it 
holds for the general case of a diagonal Hamiltonian. 



III. TREMBLING CAVITIES 

Now we are in the position to apply the formalism 
presented in the previous Section to a special system 
of a quantum field under the influence of dynamical 
external conditions. We consider a massless and neu- 
tral scalar field $ confined in an arbitrary and weakly 
time-dependent domain G(t) (a trembling cavity) with 
Dirichlet boundary conditions $ = at dG(t), see also 
Refs. and [l7[ . The action generating the equation 
of motion is given by 



A 



1 



dt / d 3 r(9 Al $)(^$) 



dtL. 



(29) 



G 



The Hamiltonian governing the dynamics of this sys- 
tem can be obtained in the following way: For every 
fixed time t we construct a complete and orthonormal 
set of eigen functions of the Laplacian inside the time- 
dependent domain G(t) satisfying the required Dirichlet 
boundary conditions. Owing to the time-dependence of 
the Dirichlet boundary conditions at dG(t) this operator 
becomes time-dependent as well. Hence also the proper 
eigen functions are time-dependent fi = fi(t). We as- 
sume a finite domain G(t) resulting in a purely discrete 
spectrum. The insertion of the expansion of the field 
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(30) 



into the Lagrangian in Eq. ( p9| ) reveals 

w?(i)Qj) +QiMu(t)Qj 



1 



+^QiMu{t)Mj K {t) Qk 



(31) 



Since the time-derivative of the field $ may also include 
the explicit time-derivative of the eigen functions, we ob- 
tain an additional anti-symmetric inter-mode coupling 
matrix 



M u (t)= J (Prhfj. 

G(t) 



(32) 



Furthermore, the eigenvalues of the Laplace operator 
uj(t) are time-dependent in general. By means of the 
Legendre transformation we obtain the Hamiltonian 



H = - (Pf - uj(t) Qj) + Qj Mu(t) Pj 



(33) 



where Pj and Qk are the canonical conjugated variables. 
In the following the undisturbed eigenvalues (frequen- 
cies) are denoted by (|i| T °°) — w | an d their variation 
by Aui 2 (t) = uij(t) — uj]. After the canonical quantization 
and the separation of the undisturbed part 



(34) 



the perturbation Hamiltonian may be cast into the form 



Hi(t) 



W(t) 



Q 2 j(t)ALu 2 j(t) + Qj(t)MjK(t)P K {t) . (35) 



AE{t) 
1 
2 



The term including Auij(t), i.e. AE(t), arises from the 
change of the shape of the domain G{t). This term AE(t) 
is called squeezing contribution. The inter-mode coupling 
■Mjk (t) contained in Wit) - the velocity contribution - 
arises from the motion of the boundaries. Of course, the 
undisturbed Hamiltonian can be diagonalized through a 
particle definition Hq = ujj(Nj + 1/2) employing the 
usual creation/annihilation operators 



a. j 



iP.i 



(36) 



With the aid of the equations above it is now possible to 
calculate the expectation value of the number operator 
employing the results of the previous Section. The eval- 
uation of the trace Tr{. . .} is most suitably performed 
in the basis of the Ho-eigenkets. One obtains a non- 
vanishing trace only for those terms that contain the 
same number of creation and annihilation operators for 



every mode. Inserting the particular interaction Hamil- 
tonian H^t) = AE(t) + W(t) into Eq. (|lj) generates at a 
first glance four terms. However, owing to the antisym- 
metry of the matrix governing the inter-mode coupling 
M.jj = the mixing terms vanish 



TrjiV/ 
Tr ' 



0, 



\Ni [w^t'Mo] A£(f)} =0 



(37) 



Consequently the squeezing and the velocity contribution 
decouple on the level of the quadratic response 



AN! = TrlN! 



-Tr <^ N T 



dtAE(t),g 



dtW(t),go 



dtAE(t)^ 



dt W(t) 



= ANj + ANY 



(38) 



Expressing Hi(t) = AE(t) + W(t)_ with the aid of the 
matrices iS and U as done in Eq. (fla) one observes that 
the squeezing effect manifests in the diagonal elements of 
the matrices S and U while the velocity effect generates 
off-diagonal elements only. 



A. Squeezing 

The diagonal form of the squeezing Ai?(t)-part of the 
perturbation Hamiltonian H± (t) indicates the highly res- 
onant character of this contribution. This fact simplifies 
the calculation of ANi because only one mode - the mode 
I - survives in the trace. With the abbreviation 



dt Auj(t) exp (2iuit) 



1 



(39) 



the squeezing part S s of the matrix S = S s + S v can be 



cast into the form 



(40) 



Of course, the squeezing part U of the matrix U = U s + 
U v is also strictly diagonal and hence does not contribute 
to ANi (see Eq. (|l|)) with the result that only the S s 
term is responsible for particle production. Therefore the 
squeezing contribution to the number of particles reads 



AN? 



Qf | (l + 2(JVj 
= I Qf 1 2 1 I 



exp(/3cjj) — 1 



(41) 



In accordance to the results of Sec. II H the particle pro- 
duction rate at temperature T equals the rate at zero 
temperature times the thermal distribution factor. 



G 



B. Velocity 

Due to the more complicated structure of the velocity 
term, i.e. the off-diagonal elements, the calculation of the 
number of created particles involves an additional sum- 
mation. Hence the velocity contribution may not be cast 
into such a simple form as the squeezing term. The W- 
part of the perturbation Hamiltonian can be expanded 
with the aid of the matrices 



9 y 



- I dtM JK (t) exp + uj K ]t) 




(42) 



and 



JK 




(t) exp(i[wj - uj K ]t) 



(43) 



In this case both terms, the S and the U matrices 
contribute. The U v term may even decrease the number 
of particles in a given mode, see Eq. (pT|). For all that 
the total energy increases. 



C. Moving mirror 

In order to illustrate the velocity effect we consider 
the most simple example of a single mirror in 1+1 di- 
mensions. This scenario has already been investigated 
by several authors, e.g. 1^,0, at zero temperature. In 
this case the domain G takes the form G(t) = (r/(t) ,oo) 
where T](t) denotes the time-dependent position of the 
mirror with rj(t J, —00) = r)(t f +00) = 0. The index 
I can be identified with the wave number / — > k = u)% 
which assumes all positive real numbers. As the shape of 
the domain G(t) does not change only the velocity-effect 
contributes. The inter-mode coupling matrix is given by 
(cf. 0,0) 



Mu(t) ^ M kk >(t) 



TT 



kk' 



k 2 - k' 2 



(44) 



where V denotes the principal value. The Fourier trans- 
form of a time-dependent function is denoted by a tilde: 
rj = Tr\ (see Eq. (|39|)). Using this notation the matrices 
for the evaluation of the particle number read 



Skk' 




and 




7 r](k + k') 



7 rj(k-k'). 



(45) 



(46) 



Inserting these expressions into Eq. ( Jig ) yields the fol- 
lowing results for the expectation values of the number 
operator 



AN k = 1 dk'% m + fcor (1 + ) + (n, 



H- / dfc'^^(fc~fc')| 2 ((iV fe ') o -(iV fe/o 



(47) 



Already for the most simple example the velocity contri- 
bution cannot be cast into a form being as simple as the 
squeezing term. But the formula above allows us to calcu- 
late the number of created particles within the quadratic 
response for arbitrary dynamics rj(t) and temperatures 
T. Deriving the total radiated energy from Eq. ( [47| ) 

00 

E = J dkkAN k = J dt J dt' r)(t)r)(t')ll(t-t') , (48) 


where 1Z denotes the quadratic response function, we ob- 
tain a more elucidative formula. For that purpose we 
have to perform integrations involving Bose-Einstein dis- 
tribution functions entering in (Nk)o- If we insert the 
usual expansion for those functions 



Nk) = 



exp(/3fc) 



I 00 
Ik) - 1 = E 



-nfik 



(49) 



n=l 



the wave number integration TT 1 leads to Hurwitz zeta- 
functions 

00 y 

If = I dkk m exp (ik(t - t') - n(3k) 



E 



m+l 



{ ( n 0-i[t-t']) 

ml J t-t' 
r C(m+l,l-t 



P 



(50) 



In terms of these functions the response function 1Z = 
TZ(At) = TZ(t — t') yields after the k- and ^'-integrations 



^ = 4 \+^(At) 



24i 
2i 



1 

24 



1 



+^ [~iTi5 w {At)-V 
■ZS (-^(At)+V 
-iTrS {1) {At) + V 



7T^ 

2 



At 2 

2i " 



At 3 
1 



At 2 



C(2,l- iAt/0) 
P 2 

C(3, 1 - iAt/0) 
P 3 

C(2, 1 - iAt/0) 
P 2 

C(3, 1 - iAt/0) 



P 3 



(51) 
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The five terms above correspond directly to the five terms 
in Eq. (ft7l). As one can easily check in Eq. (|4§|), only 
the symmetric part TZ sym (At) = (K(At) +K(-Ai))/2 of 
the response function IZ(At) contributes to the radiated 
energy. Symmetrizing the response function a lot of can- 
cellations occur and all divergent terms of the structure 
l/At n disappear. The resulting expression reads 

ft sym (Ai) = -L<5< 4 > (At) - -C(2)T 2 S (2 \At) , (52) 
iZn 7r 

with the Riemann zeta-function £(n) that is related to 
the Hurwitz zeta-function £(n,m) via £(n, 1) = C( Tl )- 
Rewriting this expression into the total radiated energy 
yields 

E=±-Jdtf l 2 (t) + ^T 2 Jdti 1 2 (i). (53) 

The first term describes the vacuum contribution and 
was originally obtained by Fulling and Davies || using 
the conformal invariance of the scalar field in 1+1 dimen- 
sions and has been later calculated by Ford and Vilenkin 
via a more flexible method of perturbations of Green 
functions. In both approaches the radiated energy was 
deduced by means of the point-splitting renormalization 
technique. The relevant contributions of the Green func- 
tions used in Refs. j^J^] correspond to the vacuum part 
of the response function IZ(Ai) in our derivation. 

The second term is a pure temperature effect and gen- 
eralizes the vacuum results in Refs. HQ to the density 
matrix corresponding to the canonical ensemble. The 
relation between the finite temperature correction to the 
radiated energy and the vacuum contribution is of the or- 
der 0[T 2 t 2 ] where r denotes a characteristic time scale 
of the dynamics of the mirror. 

IV. RESONANTLY VIBRATING CAVITY 

Let us now investigate the finite-temperature effects 
on the dynamical Casimir-effect in a resonantly vibrat- 
ing cavity. In order to allow for an experimental verifi- 
cation the number of motion-induced particles should be 
as large as possible. One way to achieve this goal is to 
exploit of the phenomenon of parametric resonance. It 
occurs in the case of periodically time-dependent pertur- 
bations characterized by some frequency lo. Within the 
quadratic response the number AiVj of created particles 
is proportional to the Fourier transform of the perturba- 
tion function (see Eq. j39|)). Assuming a harmonically 
oscillating disturbance the Fourier transform possesses a 
pronounced maximum at the resonance frequency. As 
a result particles with a mode frequency corresponding 
to lo will be produced predominantly. Obtaining large 
numbers may indicate that one has left the region, where 
second-order perturbation theory does apply. 



A. Rotating wave approximation 

In the case of oscillating disturbances, however, it is 
possible to evaluate the time evolution operator U in 
all orders of H\ analytically employing yet another ap- 
proximation, the so-called rotating wave approximation 
(RWA, see e.g. ||). The main consequence of the RWA 
consists in the fact that it allows for the derivation of 
a time-independent effective Hamiltonian Hf s after per- 
forming the integration over time 

J dtHi^ THf. (54) 

Let us assume that the explicit time-dependence of the 
perturbation Hamiltonian Hi(t) possesses an oscillatory 
behavior like e sin(2wi) during a sufficiently long time T, 
such that the conditions wT >1, £<1 and e loT = 0[1] 
hold. Expanding the time evolution operator U into pow- 
ers of e and loT the RWA neglects all terms of order 
0[e n (uiT) m ] if n > m holds. Since time integrations 
over oscillating functions result in smaller powers of T 
than the same integrations over time-independent quan- 
tities, within the RWA all terms including oscillations 
were omitted. Accordingly, only those terms, where the 
oscillations due to the time-dependence of the operators 
(in the interaction picture governed by Hq) and the ex- 
plicit time-dependent disturbances cancel - i.e. which are 
in resonance (n = m) - contribute in the RWA. This ap- 
proximation enables us to neglect the time-ordering T 
as well. The difference between the time-ordered and 
the original expression always contains commutators like 
[Hi(t),Hi(t')}. These quantities are always oscillating 
and therefore can be neglected within the RWA. 

B. Fundamental resonance 

For a harmonically vibrating cavity the effective 
Hamiltonian Hf can easily be calculated from the inter- 
action operator in Eq. (]3q) . Assuming harmonic time- 
dependences ~ e sin(2wi) or ~ £ cos(2wi) for both, the 
squeezing (Awf(t)) and the velocity terms (M-jKif)), 
only those terms will survive, which match the reso- 
nance conditions. These conditions are fulfilled if the 
oscillations of the operators, i.e. Qj(t) and Pxit), com- 
pensate the oscillations of the disturbances, i.e. Aw|(i) 
(squeezing) and M.jjc(t) (velocity). For the squeezing 
term the resonance condition reads luj — lu and for the 
velocity term \loj ±luk \ — 2a;, respectively. Accordingly, 
the squeezing effect creates always particles with the fre- 
quency lo provided this cavity mode does exist. We re- 
strict our further consideration to the situation, where 
the oscillation frequency lo corresponds to the lowest cav- 
ity mode, i.e. to the fundamental resonance 

lo = min{w/} = lo\ . (55) 
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The fundamental resonance frequency u>± is determined 
by the characteristic size A of the cavity, e.g. u>\ — 
V37T/A for a cubic cavity. For the lowest mode 1=1 
the resonance condition for squeezing u>i = u>\ = u) is 
satisfied automatically. Although the condition for the 
(S-term of the velocity effect \uoj + u>k\ — 2u> could be 
satisfied for J = K = 1 (we assume a non-degenerate 
ground state 1=1), this term does not contribute since 
M. j j = = 5,/,/ . Whether the resonance condition for 
the U-term of the velocity effect \ujj — ojk\ = 2w can be 
satisfied or not depends on the spectrum of the particular 
cavity under consideration. For a one-dimensional cavity 
the eigen frequencies ujj are proportional to integers and 
thus it can be satisfied leading to an additional velocity 
contribution. For most cases of higher-dimensional cav- 
ities, e.g. a cubic one, this condition cannot be fulfilled. 
Thus the velocity effect does not contribute within the 
RWA (cf. Ref. §). 



C. Squeezing operator 

In the following calculations we assume a case for which 
only the squeezing term contributes (i.e. the rather gen- 
eral case). The effective Hamiltonian can be derived im- 
mediately from the only contributing Aw^-terms 



,// // , = / dt ^M. ( tf Juzt + a je -*-rf 



ujje 





dt sin 



(56) 



Therefore the time evolution operator U coincides in the 
RWA with a squeezing operator S\ for the lowest mode 
1=1 



U = T 



exp 



csp [ — i j dtHi(t) 
f uel 



((«l) 2 -(«i) 2 ) 



Si, 



(57) 



with a squeezing parameter E = weT/2. This confirms 
the notion of the Awf-terms in the perturbation Hamil- 
tonian ( p5| ) as squeezing contribution. Having derived a 
closed expression for the approximate time-evolution op- 
erator U ~ Si we are able to calculate the expectation 
value for the number operator to all orders in Hf within 
the RWA 



N T ) w Tr {qoSINtSi} = ^ ^ 



+5n sinh 



weT 



(l + 2(A)J- (58) 



This non-perturbative result states that also at finite 
temperature the number of photons ANi created reso- 
nantly in the lowest cavity mode increases exponentially. 
The vacuum creation rate AA^ = sinh 2 (weT/2) (see 
Ref. M) gets enhanced by a thermal distribution factor. 
Since the effective Hamiltonian becomes diagonal in the 
resonance case such a behavior is consistent with the re- 



sults in Sec. II H 



D. Local quantities 

So far we have considered merely the expectation val- 
ues of global observables such as particle number and en- 
ergy. But the canonical formalism developed here is also 
capable for investigating local quantities. As an example 
we may consider the two-point function 



($(r)$(r')) =Tr{£$(r)$(r')} 



(59) 



According to the results of the previous Sections within 
the RWA the time evolution operator appears as a 
squeezing operator for the lowest mode 1=1. Expanding 
the field $(r) into the modes // yields 

($(r)$(r')) =Y, rTT {^S\QiQjSi}f I (r)f J (r'). (60) 

For I J the trace above vanishes and for J = J ^ 1 
it coincides with the undisturbed (thermal) expression. 
Hence the only interesting case is I = J = 1. In this 
situation the amplitudes Q± are squeezed by the time- 
evolution operator §i . As a result the change of the cor- 
relation function induced by the dynamics of the cavity 
can be cast into the form 



A($(r)$(r' )> = (e 23 - l) Tr {g Q Ql} 
xfi(r)fi(r') 
= (e 23 - 1) (l + 2 (iVi 



x/i(r)/i(r')/(2wi) 

x/i(r)/i(r / )/(2w 1 ), (61) 

where S again denotes the squeezing parameter. In com- 
plete analogy one obtains the change of the correlation 
of the field momenta 



AOl(r)n(r')) = (e~ 2a - 1) 1 



exp(/3a>i) — 1 / 
xfi(r)fi(r')wi/2. (62) 

Note, that the canonical variables Qi and the momenta 
Pi transform in an opposite way under squeezing: Qi — > 
e~Qi whereas Pi — > e~ a Pi. 
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These ingredients enable us to calculate the expec- 
tation values of the energy-momentum tensor T^v = 
d^&dvQ — g^d p $d p &/2. The expression for the change 
of the energy density reads 

a ™4( 1 + ^ft) 

xP-lj^fM + le- 25 -!)/^] . 

(63) 

The stress tensor Ty consisting of purely spatial com- 
ponents of T^y can be calculated in a completely ana- 
logue manner - one would obtain additional terms like 
difidjfi. This quantity can be used to deduce the me- 
chanical properties, e.g. the pressure, of the radiation 
field inside the cavity. The expectation value of the en- 
ergy flux density Tbi vanishes within the RWA, since the 
change of the energy is always one power of the vibration 
time T lower than the energy itself. 

The above expression can be used to deduce the spa- 
tial distribution of the energy created by the dynami- 
cal perturbation of the cavity. Since the squeezing pa- 
rameter S is proportional to the vibration time T, the 
first term at the r.h.s. dominates for large time dura- 
tions T. In this situation the produced energy density 
behaves as (V/i) 2 . For Dirichlet boundary conditions the 
eigen functions vanish at the boundary but their deriva- 
tives usually reach their maximum value there. In the 
center of the cavity the lowest eigen function assumes its 
maximum and - consequently - its derivative vanishes. 
Ergo the energy density is concentrated near the bound- 
aries of the cavity in the case under consideration. 

However, this assertion crucially depends on the im- 
posed boundary conditions. For Neumann conditions the 
behavior is actually opposite. In order to analyze the sit- 
uation for the more realistic photon field one has to con- 
sider the corresponding eigen modes for this case. The 
modes of the electromagnetic field are more complicated 
and cannot be reduced to purely Dirichlet or Neumann 
conditions in general. As demonstrated in Ref. |28|| , one 
may divide the electromagnetic eigen functions into TE 
and TM modes. E.g., for a cylinder with perfectly con- 
ducting walls the energy in the lowest TM mode is con- 
centrated near the shell and vanishes at the symmetry 
axis, see p8|| . In addition, the eigen frequency of this 
particular mode does not depend on the separation of 
the bottom and the top surface. Hence it cannot be ex- 
cited resonantly by a vibration of the piston but only by 
a change of the radius. In contrast, the lowest TE mode 
does depend on the separation and thus may be excited 
by a vibration of the piston. In summary, the conclu- 
sions concerning the energy distribution of the photon 
field crucially depend on the particular mode that is ex- 
cited by the cavity vibration. 



E. Discussion 

In order to indicate the experimental relevance of the 
calculations above one may specify the characteristic pa- 
rameters. For room temperature ~ 290 K, which corre- 
sponds to thermal wave lengths of about 50 /xm and con- 
sidering a cavity of a typical size A « 1 cm one obtains a 
thermal factor (1 + 2(Ni )o) = O[10 3 ]. As a consequence 
the number of photons AA^i created by the dynamical 
Casimir effect (after the vibration time T) at the given 
temperature T will be several orders of magnitude larger 
in comparison with the pure vacuum effect at T = 0, see 

S- . 

This enhancement occurring at finite temperatures 
could be exploited in experiments to verify the phe- 
nomenon of quantum radiation as long as back-reaction 
processes (and losses, etc.) can be neglected. Of course, 
one has also to take into account the number of photons 
(iVj )o present at the temperature T and their thermal 
variance 




The latter quantity reflects the statistical uncertainty 
when measuring the number of photons at a given tem- 
perature. In order to obtain a number of created particles 
AA^i which is not much smaller than the corresponding 
thermal variance y/ a\ (N\ ) one has to ensure conditions 
that will lead to a significant vacuum effect as well. This 
implies that the argument of the hyperbolic sine func- 
tion in Eq. (|5§|), i.e. the squeezing parameter 5 = weT/2 
should be at least of order one. An estimate of the max- 
imum value of the of the dimensionless amplitude of the 
resonance wall vibration e max < 10~ 8 is given in Ref. [||. 
For a characteristic size of the cavity of about 1 cm cor- 
responding to a fundamental frequency u> ~ 150 GHz the 
squeezing parameter EE approaches 1 after several mil- 
liseconds. It still remains as a challenge whether or not 
the requirement 5 = 0[1] could be achieved in a realistic 
experiment. But - provided an experimental device for 
generating a considerable vacuum contribution becomes 
feasible - there will be a strong enhancement of the dy- 
namical Casimir effect at finite temperatures. 

V. DIELECTRIC MEDIA 

The previous Sections were devoted to the investiga- 
tion of a scalar field confined in a trembling cavity with 
Dirichlet boundary conditions simulating perfect conduc- 
tors. Although the effect of quantum radiation has not 
yet been verified conclusively in an according experiment, 
a resonantly vibrating cavity is expected to provide one 
of the most promising scenarios for this aim. Of course, 
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the assumption of perfectly conducting walls of the cavity 
is an idealization. One possible step towards a more re- 
alistic description is to consider a dielectric medium with 
a finite permittivity e. Of course, one may also take into 
account the permeability, see e.g. [^?|. But for reasons 
of simplicity we restrict our further considerations to a 
purely dielectric medium. 

The quantum radiation generated by a moving body 
with a finite refractive index was studied in Ref . [[30| , for 
example. More generally, one may consider a medium 
with an arbitrary changing permittivity e(t, r) and a lo- 
cal velocity field v(t,r). Again these properties of the 
medium are treated classically, i.e. as an external back- 
ground field. As the quantum field propagating in this 
background we consider the electromagnetic field. For 
non-relativistic velocities of the medium the Lagrangian 
density governing the dynamics of the electromagnetic 
field is given by (see e.g. 29 3m) 



C = i (eE 2 - B 2 ) + (e - 1) v ■ (E x B) 



(65) 



The particle definition for this vector field requires ad- 
ditional considerations. Since it is described by a gauge 
invariant theory it possesses primary and secondary con- 
straints, see e.g. In Ref. |^9| these gauge problems 
are solved by virtue of the reduction of variables. How- 
ever, other procedures, e.g. the Dirac quantization, lead 
to the same results Q . 

There have been various efforts to discover effects of 
quantum radiation for such dynamical dielectrics: One 
interesting idea goes back to Schwinger Q, who sug- 
gested to explain the phenomenon of sonoluminescence 
by this mechanism. Sonoluminescence means the con- 
version of sound into light. In an according experiment 
one generates sound waves in a liquid (water) in such 
a way that tiny oscillating bubbles emerge. Under ap- 
propriate conditions these bubbles emit light pulses, see 
e.g. j3j| and references therein. In spite of the consider- 
able amount of work and the controversial discussions in 
order to clarify the relevance of quantum radiation with 
respect to sonoluminescence, see e.g. pij , ]3q] , |37|] and 
also p9| , p0| , there are still open questions, since the dy- 
namics of the bubble and the behavior in its interior are 
not known sufficiently. We shall return to this point later 
on. 

In view of the Lagrangian density in Eq. ( |65| ) the dy- 
namical properties of the medium enter in two terms: In 
analogy to the cavity example we may distinguish be- 
tween the squeezing effect due to a varying permittivity 
e(t, r) and the velocity effect governed by v(t 7 r). In the 
following we consider situations where the squeezing term 
gives the dominant contribution (see Ref. [£9| , Section V) 
and neglect the velocity field (v = ). In contrast to the 
cavity example the squeezing term of the perturbation 
Hamiltonian will not be diagonal in general 



Hi(t) 



d 6 r- 



1 



1 



1 



2 \e(t,r) 



n (t,r) 



= / <FV0(f,r)II (t,r). 



(66) 



II = E denotes the canonical momentum density asso- 
ciated to the vector potential A, see e.g. |^] and flT7| . 
9(t, r) symbolizes the deviation of the permittivity e(t, r) 
from its asymptotic value eoo = e(\t\ t oo, r) — e(t, \r\ f 
oo). 

The diagonalization of the undisturbed Hamiltonian 
via a particle definition can be achieved with photons 
labeled by / = {kv} where k denotes the wave vector of 
the photon and v counts the two possible polarizations. 
Within this basis the S and U matrices assume the form 

&kv,k'v' — ^ yS-kv ■ ek'v') 

x / d 4 x9(x) exp {i(k + k')x) 



{e-kv ■ e k 'u') 6{k + ti) 



where 13 denotes the quantization volume, and 



y/UJk^k' 





(e-kv -e fe v') 0(k-k'), 



(67) 



(68) 



respectively. In complete analogy to the cavity we may 
calculate of the quadratic response of the number opera- 
tor 

AN k „ = l<w,'| 2 (i + (^w) + (#*") ) 

h'v' 

I^W-'l 2 ((^W) o - (#*«,) ) • (69) 



k' 

It is again possible to recognize the thermal corrections 

I . In Ref. 129] we gave 



to the vacuum effect ^ \Sk v ,k'v' 

k'i/' 

a general proof that for massless and non-self-interacting 
bosonic fields at zero temperature the spectral energy 
density e(w) created by smooth and localized distur- 
bances behaves as e(uj) ~ uj a for small u>. As one can 
easily check in the equation above this is no longer valid 
in general at finite temperatures due to the Boltzmann 
distribution function that becomes singular with 1 jio for 
small lo. 



A. Small R expansion 

The structure of Eq. (^9|) is too complicated for a 
general discussion of the physical properties of the in- 
duced quantum radiation by means of simple expressions. 
For that purpose it is necessary to use some approxima- 
tions. One possibility is to assume that the region where 
the permittivity e(t,r) differs from its asymptotic value 
£oo = e(|^| T oo, r) = e(<, \r\ f oo) is very small. This 
assumption can be used to expand the Fourier transform 
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of the perturbation function 9(t, r) in powers of i?, where 
R denotes a characteristic length scale of the disturbance 



8(k) 



dt e" 



r e 



k-r 



0(t,r) 



o / dte lu)t V(t) + 0[R 4 



6 V(lj) + O[R 4 



(70) 



For the lowest (volume V ~ R 3 ) term of this expan- 
sion it is possible to calculate the associated radiated 
energy in close analogy to the moving mirror example. 
But in the case under consideration the evaluations have 
to be accomplished in 3+1 dimensions which leads to 
additional scale factors k 2 and k' 2 due to the d 3 k- and 
d 3 /c'-integrations. This results in the occurrence of Hur- 
witz zeta-functions of higher order £(4, 1 — i(t — t')j (3) 
and C(5, 1 — i(t — t')//3) and therefore in higher powers 
of the temperature T. After some calculations the lowest 
order terms in R and T of the total radiated energy yield 



E 



2^ 



3 

dt V 



5 • 7 

2 



C(4)^T 4 I <HV 2 



(71) 



In analogy to the moving mirror example in Sec. Ill C 



the first term describes the (lowest) vacuum contribution, 
which was already obtained in Ref. 29 1 , whereas the sec- 
ond term represents the (lowest) temperature correction. 
But in contrast to the moving mirror in 1+1 dimensions 
the lowest thermal correction increases with T A in this 
scenario owing to the additional k, ^'-integrations. 

If we would have the exact data for the oscillating bub- 
ble we were able to evaluate the number of photons gener- 
ated by the quantum radiation and so we could quantify 
the contribution of this mechanism to the phenomenon of 
sonoluminescence - under the assumptions made. But as 
these data are not known yet with a sufficient accuracy, 
this question remains unsolved at this stage. 



B. Large R limit 

Now we consider the opposite situation and assume 
that the permittivity changes over very large volumes in 
the same way. In such a scenario the disturbance function 
9 becomes nearly position-independent 6{t,r) ~ 9{t). In 
this limit it is also possible to simplify Eq. (p9) since the 
<i 3 r-integrations in Eqs. ( |67| ) and ( p8| ) produce 8 3 (k±k')- 
distributions and therefore the mode integrations break 
down. As a result the expression for the quadratic re- 
sponse of the number operator obeys a structure similar 
to the squeezing term in the cavity example 



AJVfc„ = \S kUt{ _ k> \ [l + 2(N, 



(72) 



To establish the analogy once more we note that a per- 
turbation like 8(t, r) — e sin(2o;i) generates also a gener- 
alized squeezing operator 



17(H) = exp 



E 

v . I fc I — [ 



(P 4 f 



h.c. 



(73) 



similar to the resonantly vibrating cavity, see also 
pSfl . However, there is a crucial difference between the 
medium and the cavity: In a closed cavity there exist only 
particles with special discrete frequencies (eigenvalues 
of the Laplace operator). In contrast, for the dielec- 
tric medium without boundary conditions all positive 
frequencies are occupied by photons. Hence one has to 
vibrate a (finite) cavity with a special (resonance) fre- 
quency in order to create particles while in a medium 
the frequency may be arbitrary. 

Investigating the two-photon correlation Ckvk'v' a t 
zero temperature one observes that in this case 
the photons are most probably emitted back-to-back: 
Ckv,k'v'{T = 0) ~ 8{k + kJSw'. At finite temperatures 
the second term in Eq. (EQ) gives raise to an additional 
negative correlation which is isotropic, i.e. it does not de- 
pend on the directions of propagation of the two photons. 
However, the anisotropic and temperature-independent 
back-to-back correlation represents one possibility to dis- 
tinguish between the photons arising from the quantum 
radiation and the purely thermal radiation. This observa- 
tion (sec |3£j ) might perhaps be used to clarify the origin 
of the photons (i.e. the underlying mechanism) within 
the phenomenon of sonoluminescence. 



VI. FRIEDMANN-ROBERTSON- WALKER 
METRIC 

In the previous Sections we focused our attention on 
mirrors represented by Dirichlet boundary conditions and 
on dielectric media. Now we are going to apply the 
canonical formalism developed there to yet another sce- 
nario - where the gravitational field generates quantum 
radiation. 

According to the commonly suggested scenario the cos- 
mological evolution starts at a stage of high tempera- 
tures. It is generally believed that the back-reaction 
of the cosmological particle production onto the gravi- 
tational sector yields a significant contribution. Conse- 
quently, it will be important to calculate the tempera- 
ture effects that could affect the cosmological dynamics 
at very early stages. 

Let us consider the minimally coupled massless scalar 
field propagating in the conformally flat Friedmann- 
Robertson- Walker space-time, see e.g. Ref. |39| for a re- 
lated calculation at zero temperature. 

The Friedmann-Robertson- Walker metric represents a 
solution of Einstein's equations for a homogeneous and 
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isotropic distribution of matter and describes an expand- 
ing (or contracting) universe. Depending on the density g 
of the matter (if we omit the cosmological constant) there 
exist three different branches of this solution: For den- 
sities g exceeding the critical density g c one obtains the 
closed elliptical universe, which eventually re-collapses. 
For g < g c one is lead to the open hyperbolic universe and 
for g = g c we get the open conformally flat Friedmann- 
Robertson- Walker space-time. In contrast to the first 
case (g > g c ) the other scenarios (g < g c ) imply an eter- 
nally expanding universe. 

In order to specify the correct value of the density g 
one has to deal with the problem of the unknown dark 
matter. In view of the present status of the observations 
it might well be possible that the density g is indeed 
close or equal to the critical value g c , which is connected 
to the Hubble constant. In any case, for small space- 
time domains the conformally flat Friedmann-Robertson- 
Walker space-time should be a good approximation in 
any case. In terms of the conformal coordinates (t, r) 
the corresponding metric is given by 



ds 1 



Q 2 (t) (dt 2 -dr 2 ) 



(74) 



where fi 2 (i) denotes the scale factor governing the Hub- 
ble expansion. It describes the change of the measure of 
length and time scales during the cosmological evolution, 
e.g. inducing the cosmological red-shift. 

However, the following calculations will become easier 
if we introduce a slightly different time coordinate t —> t 
with 



dr = n~ 2 dt . 



(75) 



see also Ref. [|39|. In terms of the new time coordinate 
the metric can be cast into the form 



ds 2 



n b (r)dT 2 - n 2 (r)dr 2 



(76) 



The action for a minimally coupled massless scalar field 
propagating in this particular curved space-time reads 



A 



d 4 x- 



gd^g^d^ 
i / ' dr J d 3 r ($ 2 - fi 4 (V$) ; 



(77) 



As the advantage of the new time coordinate r we observe 
the cancellation of the scale factor in front of the <I> 2 -term. 
Consequently the equation of motion assumes the simple 
form 



<9 2 $ 

W 1 



fi 4 (r)V 2 $. 



(78) 



After the usual canonical quantization procedure the 
Hamiltonian can be cast into the form 



H(r) 



1 



d 3 r(n 2 + f2 4 (r)(V$)' 



(79) 



One ob serve s a close similarity to the large R limit in 
Section VB . As we shall see below, this similarity holds 
also for the number of created particles. 

In complete analogy to Sec. HI it is possible to diag- 
onalizc this time-dependent Hamiltonian by an expan- 
sion of the field $ into a complete set of orthogonal 
eigen functions of the Laplace operator. Owing to the 
spatial homogeneity and isotropy of the space-time there 
is an ambiguity concerning the selection of such a basis 
set. Here we choose the eigen functions to be completely 
time-independent, fj = fi(r). As a consequence, adopt- 
ing this expansion 



(80) 



the resulting modes Qi(t) do not obey any inter- mode 
interaction due to the spatial integration and the orthog- 
onality and time-independence of the eigen functions //. 
Hence the time-dependent Hamiltonian is diagonal 



(81) 



where —to 2 denote the time-independent eigenvalues of 
the Laplacian. Now we may use the outcome of Section 
[I H and we arrive at 



AJVJ = AN, 



T=0 



1 



exp(/3w) 



1 



(82) 



It should be mentioned here that the particle number 
above is - strictly speaking - merely a formal quantity 
since it does not describe particles in a well-defined and 
unique sense. The Friedmann-Robertson- Walker space- 
time is not time-translationally invariant and thus does 
not possess a corresponding Killing- vector. Ergo the def- 
inition of energy necessitates additional considerations. 
It is not possible to define a physical reasonable and con- 
served energy. Of course, this fact is consistent with the 
permanent particle creation. Hence the interpretation of 
the above quantity AN^ is not obvious - at zero as well 
as at finite temperatures, see also Ref. [Q . But here we 
are mainly interested in the influence of finite (initial) 
temperatures. Fortunately, the finite temperature effects 
factorize out and the problems mentioned above concern 
the pre-factor AA^J^ only. In summary we may draw 
the conclusion that - putting aside the problem of the 
interpretation of the vacuum term AiV ( J =0 - the particle 
creation in the Friedmann-Robertson- Walker space-time 
at finite (initial) temperatures gets strongly enhanced by 
a thermal factor in analogy to the resonantly vibrating 
cavity. 

It should be mentioned here that the phenomenon of 
particle creation (quantum radiation) in the conformally 
flat Friedmann-Robertson- Walker space-time can be ob- 
served merely for fields which are not conformally invari- 
ant. As counter-examples we may quote the massless 
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scalar field in 1+1 dimensions, the massless and confor- 
mally coupled scalar field in 3+1 dimensions (see also 
Ref . Q ) , and - last but not least - the electromagnetic 
field in 3+1 dimensions. Obviously the absence of any 
mass terms is essential for the conformal invariance. For 
these conformally invariant fields the equation of motion 
does not lead to any mixing of positive and negative fre- 
quency solutions within the conformally flat Friedmann- 
Robertson- Walker metric and thus no particles are cre- 
ated. 



VII. SUMMARY 

Calculating the number of particles created by dynam- 
ical external conditions we found that for the case of 
a completely diagonal Hamiltonian the number of pro- 
duced particle at finite temperature equals the analogue 
quantity at zero temperature times a thermal factor. 

Focusing on the scenario of a resonantly vibrating cav- 
ity we were able to derive within the RWA the effective 
perturbation Hamiltonian which turned out to be diag- 
onal. As a consequence we observe an enhancement of 
the dynamical Casimir effect at finite temperatures as 
described above. 

In contrast to this non-perturbative result the finite 
temperature corrections to the energy radiated by a sin- 
gle moving mirror in 1+1 dimensions was calculated 
within response theory. 

In close analogy we derived the energy of the photons 
generated by an bubble with an oscillating radius within 
a dielectric medium. 

Finally we investigated the particle production within 
the Friedmann-Robertson- Walker universe at finite tem- 
peratures - where the Hamiltonian can again cast into a 
diagonal form. 



VIII. CONCLUSION 

As a main result of this article we have presented a the- 
oretical description of quantum radiation at finite tem- 
peratures by generalizing the canonical approach devel- 
oped earlier 0JT^,|2| : |l|Q . The major advantage of 
this formalism is its generality and flexibility. 

Depending on the characteristic scales associated to 
the perturbation the effects of finite temperatures repre- 
sent potentially significant contributions to the quantum 
radiation and hence should be taken into account for re- 
alistic estimations of this striking effect. 

This observation may be interpreted in the following 
way: Not only the vacuum fluctuations but also the ther- 
mal excitations arc converted into real particles by the 
influence of the dynamical external conditions. 



IX. DISCUSSION 

For the perhaps most interesting scenario in view of 
a possible experimental verification of the dynamical 
Casimir effect - the resonantly vibratin g cav ity - we spec- 
ified some relevant parameters in Sec. IV E . 

But as it became evident in Section |v|, a cavity filled 
with a dielectric medium with an resonantly oscillating 
permittivity e(t) = + Ae sin(2wi) generates quite sim- 
ilar effects. Identifying the dimcnsionless amplitude e of 
the vibration of the cavity in the first case with the rel- 
ative change of the permittivity Ae in the second case 
the set of relevant parameters is completely equivalent in 
both situations. 

It turns out that the thermal factor enhancing the dy- 
namical Casimir effect is of order 10 3 . But this enormous 
amplification should be contrasted to the thermal vari- 
ance of the number of particles present initially which 
may complicate the measurement of the number of pro- 
duced particles. For the relevant temperature regions 
both terms are of the same order of magnitude. 

As a consequence the finite temperature effects do not 
necessarily generate difficulties concerning the experi- 
mental verification of the dynamical Casimir effect. From 
our point of view there is no need to perform an experi- 
ment at low temperatures - which might be much more 
involved than one at room temperature. 



X. OUTLOOK 

There are several possibilities to extend the presented 
calculations to more general scenarios: Firstly, one may 
drop the idealization of a perfectly conduction cavity and 
take losses into account p3|j . In analogy, it would be in- 
teresting to study the effects of a non-trivial dispersion 
relation of a dielectric medium in view of the phenomenon 
of quantum radiation. This might be especially interest- 
ing for the investigation of the dielectric black hole ana- 
logues, see e.g. Q The investigation of the phenomenon 
of quantum radiation for interacting fields obeying non- 
linear equations of motion is rather challenging. 
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